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^ ■ To investigate a recent proposal that J/ip production in ultra-relativistic nuclear collisions 
0> . is of thermal origin we have reanalyzed the data from the NA38/50 collaboration within 
^ [ a thermal model including charm. Comparison of the calculated with measured yields 
I demonstrates the non-thermal origin of hidden charm production at SPS energy. However, 
^ ■ the ratio tp' /{J/ip) exhibits, in central nucleus-nucleus collisions, thermal features which 
O ' lead us to a new interpretation of open charm and charmonium production at SPS energy. 
Implications for RHIC and LHC energy measurements will be discussed. 

C ■ 

The suppression of J/ip mesons (compared to what is expected from hard scattering 
models) was early on predicted |I| to be a signature for color deconfinement. Data for S- 
^ I induced collisions exhibited a significant suppression but systematic studies soon revealed 
c3 I that such suppression exists already in p-nucleus collisions and is due to the absorption in 
(normal) nuclear matter of a pre-resonant state consisting, e.g., of a color singlet ccg state 
that is formed on the way towards J/ip production. The situation has been summarized 

The newest data for Pb-|-Pb collisions now exhibit clear evidence for anomalous absorp- 
tion beyond the standard nuclear absorption expected for such systems. The most recent 
results are summarized in The observed anomalous suppression is not explained in 

conventional models where the charmonia are broken up by interactions with co-movers 
as discussed in |0. For a discussion of the present status of J/ip suppression and its 
understanding in terms of phenomenological models see 

However, it was recently conjectured [§,|l3] that J/ip production is of thermal origin 



and exhibits no direct connection to color deconfinement. Since the charmonia are heavy 
mesons with masses much larger than any conceivable temperature, thermal production 



would be a big surprize. On the other hand, substantial evidence now exists |TT|--[17| that 
hadron production (other than charmonia) in ultra-relativistic nuclear collisions proceeds 
through a state of chemical equilibrium near or at the phase boundary between hadron 
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matter and quark-gluon plasma. For a review of the implications for quark-matter physics 



see 



To shed more light on the situation we have modified the thermal model used in [jTT| to 
include charmed hadrons. We will first present an overview of the relevant modifications 
along with a comparison of resulting thermal yields for ^/ip mesons and of the if) /{J/tp) 
ratio with experimental results from the NA38/50 collaboration. These results and a 
comparison to charm production in hard scattering models suggest a new approach to- 
wards understanding charmonium production in ultra-relativistic nuclear collisions, which 
is discussed in the following. The implications for collider experiments at RHIC and LHC 
will also be summarized. 



The present statistical model is based on the use of a grand canonical ensemble to 
describe the partition function and hence the density of the particles of species i in an 
equilibrated fireball: 

with Ui = particle density, Qi = spin degeneracy, h = c = 1, p = momentum, and E = 
total energy. The chemical potential including charm degrees of freedom is written as 
f^i = fJ'sBi + fJ^sSi + fJ'i^If + fJ'cCi- The quantities Bi, Si, If, and Cj are the baryon, 
strangeness, three-component of isospin, and charm quantum numbers of the particle of 
species i. The temperature T and the baryochemical potential fiB are the two independent 
parameters of the model, while the strangeness chemical potential ns, the charm chemical 
potential fic, and the isospin chemical potential fij^ are fixed by strangeness, charm, and 
charge conservation. In addition, the volume V of the fireball is determined by baryon 
conservation via the relation UbaryorN = ^part, where Np^rt denotes the number of nucleons 
participating in the collision and Ubaryon is the net baryon density computed in the thermal 
model. 

In addition to the standard hadronic mass spectrum of 191 hadrons as used in we 
have added mesons and baryons with open and hidden charm. Specifically, open charm 
particles included are: D''", D^, D°, D^, Sc, A*, A^ S^, A*, along with the charmonia 
Vc J/^, Xo, Xi, X2, 

Since the inclusion of these hadrons will modify the rest of the hadron yields only at 
the sub-percent level, we will use, for the following investigations, the temperature T= 
168 MeV and baryon chemical potential fiB = 266 MeV as established for central Pb-|-Pb 
collisions at SPS energy in [^. This leads to fis = 71 MeV, fij^ = - 5 MeV, (both as in 
[p]]), and /ic = - 65 MeV. Using a volume of 3085 fm^ determined by baryon conservation 
for central Pb+Pb collisions corresponding to 400 participants, we have compared, in 
Fig. m the predictions of this thermal model with the data from NA50 |21,22|, as recently 



analyzed by J. Gosset et al. |2^ . Here the data for J/ip multiplicities are plotted vs Impart- 
For this conversion we used the relation between transverse energy and impact parameter 
as given for the 1995 NA50 data in |2^J2^ and the connection between impact parameter 



and Npart as established in |2J] 



The predictions of the thermal model are represented by the dashed line, where we have 
made use of the fact that, within the range of applicability of the thermal model, all yields 
scale proportional to the volume, i.e. proportional to the number of participants. We note 
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that, even for the most central colhsions (where i /ifj production should be most suppressed 
the measured i/ip yield is underpredicted by the thermal model calculations by 
more than a factor of 2 and the discrepancy is about a factor of 3 for Npa^t = 200. To 
compensate this factor of 3 discrepancy by an increase in temperature^ would require a 
temperature of T = 180 MeV, a value not compatible with that determined from hadron 
production yields [|T^], where a temperature range of 168 ±4 MeV (including systematic 
uncertainties) was established. 

We further note that it is highly doubtful that full chemical equilibration can be reached 
for charmed hadrons at SPS energies, either in a hadronic or a quark-gluon plasma sce- 
nario. Cross sections for open charm production among hadrons are in the sub-/ib level 
for relevant (thermal) energies, as is estimated from [^]. At such cross section levels the 
equilibration times for charm should exceed those for strangeness, where production cross 
sections exceed 100 /xb, by more than 2 orders of magnitude. Taking into account that 
strangeness equilibration times in a hadronic fireball exceed 50 fm/c p6|, chemical equi- 
librium for charm in the hadronic sector can be ruled out. Even in a quark-gluon plasma, 
where cross sections for charm production are much larger, thermal production is small. 
Assuming a charm quark mass of 1.5 GeV and an initial temperature of 300 MeV, very 
high for SPS energies, Redlich has estimated in a parton cascade approach p7|j28[] that, 
at hadronization (with Tc = 160 MeV), the number of thermal cc pairs is less than 0.01 
in central Pb+Pb collisions, lower by a factor of 40 of what would be needed to explain 
the data. Furthermore, at T = 160 MeV, before the start of the mixed phase, we estimate 
the volume of the plasma phase to be Ypiasma = t^R^t = 950 fm^, assuming as in |]27 



a lifetime of the plasma phase of 6.7 fm. The number of charm quark pairs in chemical 
equilibrium is then N^f = 0.47, implying that in the cascade approach the parton gas is 
undersaturated in charm by about a factor of 50! 

In Table |I] we present a summary of the results from the thermal model calculation 
for mesons and baryons with open and hidden charm^. It is interesting to compare these 
numbers with predictions for the production of hadrons with open charm in hard collisions. 
Results obtained by PYTHIA calculations following ||25| are shown in Table |^. We note 
that, somewhat surprizingly, the yield of directly (via hard collisions) produced charm 
is slightly larger than that produced if charm is in full chemical equilibrium at T= 168 
MeV. This implies that one cannot, under any circumstances, neglect direct production. 
On the other hand, ratios of charmed meson yields differ significantly from the thermal 
to the direct scenario, as is obvious from Table |I] and Table |^. 

Another interesting point concerns the ip /{J/ip) ratio. As is well known |29|, this ratio 
is, in hadron-proton and p-nucleus collisions, close to 12 %, independent of collision sys- 
tem, energy, transverse momentum etc. In the thermal model, the ratio is 3.7 %, including 
feeding of the J/ip from heavier charmonium states. A temperature of about 280 MeV 
would be necessary to explain the ratio found in pp and p-nucleus collisions in a thermal 
approach. Clearly, J/ip and production in pp and p-nucleus collisions are manifestly 



^From T = 168 MeV to T=178 MeV, the 3/ip yield increases by a factor 2.7, the (J/V')/7r~ yield increases 
by a factor of 2, and the ip /{3/4>) ratio by a factor 1.2. 

•^The thermal fluctuations about these mean values n*'''^''™ are Poisson distributed. This implies that per 
collision the variance equals n*'"^'"™. For 10^ Pb+Pb collisions, the thermal prediction is, consequently, 
that 200 ±14 3/4> mesons will be produced. 
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non-thermal. This was previously realized by Gerschel Similar considerations apply 
for the X states. In fact, feeding from xi to J/'?/' is less than 3 % if the production ratios 
are thermal. 

The evolution with participant number of the ip /{J/ip) ratio in nucleus- nucleus colli- 
sions is presented in Fig. |^. The data are from the NA38/50 collaboration [^l|-|3^. With 
increasing Np^rt the ip /{J/ip) ratio drops first rapidly (away from the value in pp col- 
lisions) but seems to saturate for high Npart values at a level very close to the thermal 
model prediction, both for S+U and Pb+Pb collisions 

This surprizing fact along with the previous observations leads us to propose a new 
scenario for J/ip production. We assume that all cc pairs are produced in direct, hard 
collisions, i.e. in line with previous considerations we neglect thermal production. For a 
description of the hadronization of the c and c quarks, i.e. for the determination of the 
relative yields of charmonia, and charmed mesons and baryons, we employ the statistical 
model, with parameters as determined by the analysis of all other hadron yields |[Tl|]. The 



picture we have in mind is that all hadrons form within a narrow time range at or close 
to the phase boundary. 

Since the number of directly produced charm quarks deviates from the value determined 
by chemical equilibration, we introduce a charm enhancement factor Qc by the requirement 
of charm conservation. This leads to: 



Nf^^ = lgcV{j2 <™ + + 9iv{j2 + ■■■■ (2) 

i i 

Via this equation the thermal yields (thermal densities times volume V) are adjusted 
to the yield of directly produced charm quark pairs. The remaining terms in eq. (2) 
of second order and third order in are completely negligible. Of course this equation 
makes only sense as long as N^^'"^'^* is much less than the number of up, down, and strange 
quarks in the fireball, a relation which is well fulfilled up to the highest (LHC) energies 
considered here. 

The number of D and J ftp mesons are then enhanced relative to the thermal model 
prediction by factors gc and g^, i.e. 



Nd = gcVnT"^ and Nj/^ = g.^Vn*);™. (3) 

Using this approach of direct production and statistical hadronization we have recalcu- 
lated the yield for J/?/) mesons. For Np^ri = 400 in Pb+Pb collisions the directly produced 
number of charm quark pairs is N'^^''^'^^ = 0.173 per collision, using the PYTHIA parame- 
ters of [^]. This implies gc = 1.38. The resulting J/ip yield per participant is plotted, in 



Fig. 1^ as a function of Impart Q- 

Without introduction of new parameters this approach describes the measured yield 
for mesons very well for those collision centralities where also the ip /{J/ip) ratio 

^For the ^part dependence of open charm we use the form deduced from the dependence on transverse 
energy of the Drell-Yan yield, as determined by p3| . Somewhat surprizingly, this scales approximately 
linearly in Npo^t, not like N^^Jj as expected for hard scattering. 
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is well described (see Fig. H). The interpretation of these results is as follows: at SPS 
energies charm is produced directly in nuclear collisions at the rate expected by straight 



extrapolation from what is known about charm in pp collisions [25|. From these charm 



quarks all charmed hadrons are formed by statistical hadronization. 

We would like to point out, however, that the present approach is rather schematic. 
First, the absolute yield of charm in nucleus-nucleus collisions is not well known. From 



25| we deduce an uncertainty of the K factor of about 2. Furthermore, open charm pro- 



duction might be enhanced by a factor of 3 in Pb+Pb collisions over expectations from 



pp collisions, as has been conjectured in |24]. Finally, inclusion of additional charmed 



mesons (in particular D*'s), will modify the absolute yield of thermal charm. A quantita- 
tive description will require a measurement of the open charm yield in Pb-|-Pb collisions. 

An interesting question was raised after the first submission of this paper, namely 
whether the canonical or grand canonical partition function should be used In the 
present approach, we have chosen the simpler grand-canonical ensemble. We remark that 
the Impart dependence of the J/ip yield may help to decide among the different approaches. 

In the present approach all J/tp mesons result, for the most central collisions, from 
the statistical hadronization of the directly produced charm quarks. Directly produced 
J/ijj mesons are (i) not formed before the reaction proceeds into a plasma phase or (ii) 
effectively destroyed during the plasma phase by, e.g., a color screening mechanism as 
proposed in In either case the current interpretation requires the existence of a 
deconfined phase during the collision. We remark here that Kabana has recently 
argued that coalescence of charm quarks is the source of J/ip mesons in nuclear collisions. 
While in spirit this is similar to our approach, this scenario assumes an enhancement of 
open charm which increases with ^part, very different from the present conclusions. In the 
context of the above arguments it would be very important to get a direct measurement 
of open charm production in nucleus-nucleus collisions. 

If the present scenario is correct, the consequences for quarkonium production at collider 
energies could be significant. For RHIC energies near central rapidity, e.g., the number 
of directly produced charm quark pairs per unit rapidity is |^ dN'^^'^'^^ / dy =1.0, leading 
to gc = 6.8 and dN^jfl^^/dj = 10~^, close to the unsuppressed value expected from hard 



collisions If the directly produced charm quarks hadronize statistically, as is implied 
at SPS energy, we would predict no J/ip suppression at all at RHIC energies, even though 
there are no J/ip mesons during the plasma phase. At LHC energy, the current approach 
would actually predict a significant enhancement of charmonia over the value expected 
for direct production. Of course, the underlying assumption is that the momenta of 
the charm quarks are close to thermal near the critical temperature, i.e. thermal (but 
not chemical) equilibration is required. Present SPS data are not at variance with such a 
scenario, since the measured transverse momentum spectra for J/ip mesons [§^,0 exhibit 
thermal shapes with inverse slope constants around 230 MeV, as expected for a heavy 
particle which participates little in the transverse flow build-up during hadronic expansion. 
Whether this thermalization will also take place at collider energies is an interesting open 
question. 



We finally note the difference of the present approach with that described in [pS 
where secondary charmonium production during the mixed phase is calculated under the 
assumption that all charm quark pairs end up in D mesons after hadronization. Since, 
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even in the present approach, the yield of charmonia is small compared to the yield of D 
mesons, the calculations reported in |^8|] concerning secondary charmonium production 
should still be valid. In fact, the charmonium production yields considered there should 
be added to the present predictions. 

In summary, we have shown that assuming chemical equilibration of charm does not 
lead to a successful thermal description of available data for p-nucleus and nucleus-nucleus 
collisions. However, the experimental data at SPS energy for the ratios ip /{J/ip) and 
(J/'?/')/Npart exhibit thermal features for the most central Pb+Pb collisions. Coupled 
with the fact that direct production of charm quarks in hard scattering is close to the 
value obtained by assuming full chemical equilibrium these observations lead to a new 
interpretation of charmonium production in nucleus-nucleus collisions in terms of a direct 
production and statistical hadronization approach. This describes the SPS data well and 
suggests a possible revision of the scenario for charmonium production at collider energies. 
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Table 1 

Yields for particles with open and hidden charm at mid-rapidity for central Pb+Pb col- 
lisions as predicted by the thermal model calculation. Chemical freeze-out takes place at 
T=168 MeV in a volume of 3085 fm"' corresponding to 400 participants. The chemical 
potentials are: fis = 266 MeV, fis = 71 MeV, fic = -65 MeV, fxj., = -5 MeV. The 
effect of feeding of the J/iJj meson from heavier charmonia, including the relevant decay 
branching ratios, is shown in the last row. 



particle species 


yield/TT 


total 


D+ = DO 
D- = DO 

Ac 
Ac 


5.8-10-5 

1.3- 10-4 

6.4- 10-5 
5.8-10-6 


0.034 
0.073 

0.036 
3.3 -10-3 


J/^ 
Xi 
^' 

J/V' + x's + ^' 


3.45 10-^ 

3.5 10-s 
1.3 10-^ 

3.6 10-^ 


2.0 -10-4 
2.0-10-5 
7.6 -10-6 

2.1 -10-4 
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Table 2 

Yields for particles with open charm for central Pb+Pb collisions as predicted by PYTHIA 
calculations p5| for N-N collisions and scaled to central (Np^^ = 400) Pb+Pb collisions 
using the nuclear thickness function. 



particle species 


total yield 


D+ 


0.032 




0.036 




0.093 


DO 


0.12 


Ac 


0.034 




0.011 




0.015 
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Figure 1. Comparison of measured J/ip yield per Pb+Pb collision as a function of 
centrality with the predictions of the thermal model (dashed line). The data are from 
as analyzed in ||23|. For details see text 
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Figure 2. Comparison of the dependence of the measured ijj' / (J/V') ratio on the number 
of participating nucleons with the prediction of the thermal model. The data are from 
the NA38/50 collaboration ||3l|-|34|. See text for more details. 
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Figure 3. Comparison of the dependence of the measured {J/ip)/'^part ratio on the 
number of participating nucleons with the predictions of the thermal model (dashed line) 
and of the direct/statistical model (solid line) The data are from [^,^,^. For details 
see text. 



